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' Abstract 

O \ In this work we study two complex scalar fields coupled through a quadratic 

^ ■ interaction in 2 + 1 dimensions using the method of bilinears as suggested by 

■ Rajeevg]. The resulting theory can be formulated as a classical theory. We study 

! the linear approximation, and show that there is a possible bound state in a range 

JlLi' of coupling constants. 

> 1 Introduction 

X ■ 

;_i ' Quantum field theory has been an essential tool for the modeling of various physical phe- 
. . nomena. One of the major problems in field theory is the understanding of relativistic bound 
states. The standart way to look at field theories is via perturbation theory around the free 
theory, and bound state problems are difficult to formulate in this approach. The most com- 
mon way is to use a Bethe-Salpeter approach for the four-point function (for two particle 
bound states) and find a self-consistency condition for the bound state solution. Typically 
this requires various approximations which may break down in the highly relativistic cases. 

One of the most successful applications of this approach is within the large- iVc approx- 
imation: in his classic paper, 't Hooft obtained a bound state equation for mesons in two 
dimensions in the large-iVc where iVc refers to the color for the nonabelian SU{Nc) gauge 
theory [jl[. This leads to a singular integral equation for the possible masses of the mesonic 
excitations. This equation is expressed in terms of the wave function of the meson given as 
a function of the fractional light-cone momentum. The analysis of this integral equation in 
M shows that there are only bound states corresponding to positive eigenvalues with finite 
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multiplicity, and these eigenvalues tend to infinity. The scalar version of this model is worked 
out in 1^ using the original approach of 't Hooft and in |^ via a Hamiltonian approach to 
the large-Nc limit. These relativistic equations behave in a very similar way to the standart 
't Hooft equation. In two dimensions we can generalize the Yang-Mills Lagrangian, since 
the gauge fields are not dynamical, by means of a nondynamical scalar field. The large- iVc 
limit meson bound state equation of these models have some other interesting features [Q. 
In 0, Aoki has generalized these bound state equations for bosons and fermions coupled 
via SU (Nc) gauge theory. A good presentation of many two dimensional models using the 
bilocal fields in the path integral formalism within the laxge-Nc limit is given in [0]. In this 
article several interesting bound state equations are derived, and further references are given. 

In Rajeev has formulated the large-A^^c model as a classical field theory using color 
invariant bilinears, and he has shown that the phase space of the theory is the restricted 
Grassmannian. The knowledge of the phase space allows one to make a variational ansatz 
for the baryons in this theory, which correspond to the large fluctuations of the field(see 0). 
Further details of this approach are given in |T^. Toprak and the author have extended this 
work to SO{Nc) gauge theory of bosons and fermions and obtained variants of the 't Hooft 



equation for these cases [0, |T^. The adjoint matter fields in the large- A^^^ limit yields again 
singular integral equations for possible mesonic strings, they exhibit a very similar bound 
state structure as the original model, but these equations are more complicated due to the 
fact that mesons are now color invariant strings of operators [13, [1^, T6[ . 

The two dimensional Yukawa model is analyzed within the light-cone method in W7, 



p!8| , |T9|. These models are more complicated due to nonlocal renormalization effects, it is 
possible to get an integral equation for bound states using some further approximations. 
A four dimensional extension of these ideas are given in pO]. The common feature of all 



these bound state equations is that they are singular integral equations. In the gauge theory 
cases these singular integral equations are rather restrictive in that they only allow for a 
discrete spectrum. In the other cases this is not necessarily true. There is usually a finite 
number (typically one) bound state. There are investigations in three dimensional QCD for 
the bound state equations of mesons (see the recent articlel^ll). Four dimensional realistic 
theories are very complicated since one has to deal with renormalization. The author is not 
knowledgable enough about these realistic bound state equations, some information can be 
found in the review article (see also p3| for a review of renormalization in the light-front 
point of view and some non-perturbative applications in this formalism). 

In this article we will apply a certain kind of mean field theory, which is a large-A^/ 
limit to two coupled complex bosons (we call it flavor symmetry to emphasize that it is 
not gauged). This theory is simple since it does not require coupling constant and wave 
function renormalization in the perturbation theory. Defining the scalar field around the 
free field theory may not be so interesting from a physical point of view, but we regard 
this as an interesting toy model. The linear approximation yields a bound state for the 
composite of two bosons. We will apply the methods of Rajeev |^ and formulate it as a 
classical field theory of bilinears. In this case (unlike the gauge theory case) this is only an 
approximation since the theory does not have to be restricted to this flavor invariant sector. 
To avoid repetitions we sometimes refer to our work on complex bosons in P3. There is an 
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interesting Bethe-Salpeter treatment^] of bound states in the broken phase of 0^ theory in 
in some sense this is similar to the model we work with. 



2 The model in the hght-cone and large- A^j limit 

We write down a U{Nf) invariant action for two complex scalars (pa and (pb, 

S = J d^d^<Pid,<Pa + dyld,<Pb-ml<Pi<Pa-ml4<Pb-^^^^^^ 

(1) 

We introduced a common coupling A, for the two complex fields (paa and cpba- We assume 
that the internal index a takes the values 1, ••, A^/- If we do not have the last term which 
couples the two fields the action would have a U{Nf) x U (iV/) symmetry, and the interaction 
term explicitly breaks this. The interaction may look discomforting , but it is easy to see 
that it is always positive. The classical ground state of the massive theory is when both fields 
are set to zero. This means that we can quantize the theory around its classical minimum 
by introducing creation-annihilation operators for the Fourier modes, as we will see below. 
Note that this theory in three dimensions is super-renormalizable from a perturbative point 
of view, so we do not expect any multiplicative renormalizations (when A < the theory is 
unstable, it can be defined by a perturbation theory analysis, but we restrict ourselves to 
positive values). 

Perturbatively, there is only one type of divergence after normal ordering as we will 
comment later on p6|. To apply the methods developed by Rajeev, we will use the hght- 



cone coordinates, introduce x+ = + x^) and x = -^{x'^ — x"^), and remains as the 

transverse coordinate. We choose x~^ as time (that is our evolution variable). A good review 
of light-front methods is given in , a good discussion of the scalar field in the light-front is 



also given in and in ||2P[ . The three dimensional scalar field theory has been investigated 



from different points of view in the articles |2^, |3^ . 



We use basically the same conventions in our previous work [gj, p^] . 
S = I dx-^dx-dx'{^<Pi{-2d.)d+<Pa + ^4(-2a_)a+0, - <Piiml - dl)<Pa - <Pl{ml - 9^)0, 

-^[(0l0a)' + (0l0.)'-2(0t0,)(0t^,)]). 

The action is first order in time x"*", this means that we are already in the Hamiltonian 
formalism. The Hamiltonian can be read off directly, 

H = I dx-dx' {<Pi{-dl + ml)<Pa + <pl{-dl + ml)<Pb + jUl^Paf + {cpW - 2(010,) (0^0,)]) . 

(2) 
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The quantization at equal time following Dirac gives 



[t\x-,x'),^ap{y-,y')] = -'-Spgn{x- -y-)6{x' - y'), (3) 



am grateful to E. Langmann for pointing this out to me 
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and the same rule applies for (j)iy. We recall that the field can be expanded in terms of 
creation-annihilation operators at initial light-front time (since the classical minimum is the 
zero configuration for the fields), 



2|p_ 

where we use [dp] = (To properly define everything we should assume that this expansion 
is given for (—00, — cq] U [cq, 00) at the end we take eo — ^ limit). 
The creation and annihilation operators now satisfy, 

[aaip~,Pi),a^\q-,qi)] = sgn(p_)5f - - gi], 

[baip~,Pi),b'^H'l-^(li)] = sgn(p_)5f5[p_ -gij. 

Since the fields are complex valued annihilation and creation operators are not related. We 
introduce a vacuum state |0 > for the Fock space construction, 

aa{p-,Pi)\0 >= for p_ > and a°^(p-,Pi)|0 >= for p_ < 
ba{p-,Pi)\0 >= for p_ > and 6""''(p-,Pi)|0 >= for p_ < 0. 

It is important to keep in mind that the operator aa{p-,pi) for p_ < 0, creates an antiparticle 
of momentum {—p-, —pi) and similarly for ba (which one can see by rewriting the above 
expansions in a more conventional way, by separating particle and antiparticle operators). 
We define normal ordering rules with respect to this vacuum as usual and denote it by a 
colon : (f)a---4'a we will be using for the computations the following relation, 

: a''\p-,pi)af3iq-,qi) := a""^(p_, pi)a^(g_, gi) - -^(1 - sgn(p_))(5[p_ - g_]5[pi - gi], (4) 

and exactly the same for b quanta. 

We have the following Hamiltonian in the quantized theory. 



H = j dx^dx-[:4>i{ml-dl)4>a: + :4>l{ml-dl' 

As it stands the Hamiltonian would not be a well-defined operator for finite Nf theory, we 
need to introduce mass renormalization terms which correspond in the diagramatic language 
the setting-sun diagrams]^, |3l|, |2^. When we take the large- A^j limit these counter terms 



become of smaller order, therefore the Hamiltonian as written will have a well-defined limit. 

We now define as an approximation a large- iVj limit and restrict the theory to the fiavor 
invariant sector. This is to be taken as an approximation to the full quantum theory. We 
introduce a set of fiavor invariant operators, which are directly written in the momentum 
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representation, and to simplify notation we write p to denote P-,Pi collectively, and S[p—q] 
6[p^ -q^]6[pi -gi]. 



C{p, q) 
C{p, q) 



Nf 
2 

Wf 

2 



N. 



f 



Nf 



: b"\p-,pi)ba{q-^,qi) : 

a"\p-,Pi)ba{q-,qi) 

b'"\p-,Pi)aa{q-,qi). 



(5) 



Note that C and C are just hermitian conjugates of each other. 

The idea behind the papers |TD[ is that when we take the large- A^/ limit the flavor invari- 
ant operators have smaller and smaller fluctuations, and if we compute their commutator, 
for example for A^^ with itself, we get. 



[Na{p,q),Nais,t)] 



Nf 
- (sgn(p 



(iV„(p, s)sgn{p_)S[q - r] - iV„(r, q)sgn{p_)S[p - s] 
- sgn{q_))5[p - s]5[q - r]). 



We assume that when we let N f — >■ oo there are proper large- iVj limits for these bilinears 
restricted to the flavor invariant states. As a result the theory becomes classical, the ex- 
pectation values of flavor invariant operators factorize as N f — ^ oo [|10|, |3^. Thus we may 
postulate a set of Poisson brackets for these classical variables: 



{Na{p,q),Na{s,t)} 
{Nh{p,q),Nb{s,t)} 

{C{p,q),C{s,t)} 

{Na{p,q),C{s,t)} 
{N,{p,q),C{s,t)} 
{N^{p,q),C{s,t)} 
{N,{p,q),C{s,t)} 



-2i{Na{p, s)sgn(p_)5[g - r] - Na{r, g)sgn(p_)5[j9 - s] 

— (sgn(p„) — sgn(g__))5[]9 — s\6[q — r]^ 

—2i(^Nh{p, s)sgn(p„)(5[g — r] — Nh{r, g)sgn(p„)5[p — s] 

— (sgn(p_) — sgn{q_))6[p — s]5[q — rfj 

-2t[Naip,t)sgn{q_)S[q - s] - N^is, q)sgn{t_)S[p - t] 
-(sgn(p_) - sgn(g_))5[p - t]6[q - s] 
-2isgn{q.)S[q - s]C{p,t) 
2isgn{q_)5[p - t]C{s, q) 
2isgn{p_)6[p — t]C{s, q) 
-2isgn{q_)5[q - s\C{p,t), 



and all the other Poisson brackets being zero. 

There are constraints coming from the U{Nf ) invariance, if we restrict the theory to the 



flavor invariant sector. They are very similar to the ones found in [24 



(eNa + e)2 + eCeCt = 1, 
(eNh + tf + eC^eC = 1, 
NaeC + CeNb + eC + Ce 
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and the hermitian conjugate of the last equation. These conditions are derived using similar 
techniques to Here we are using e(p, g) = sgn(p_)(5[p_— gi], l{p,q) = S[p—q], 



and a shorthand for the operator products, for example, NaC means, 

{NaC){p,s)= [dq^dqi]Na{p-,pi;q-,qi)C{q-,qi]S^,Si), (6) 



and similarly for the others. 

We will have further convergence conditions coming from the super-renormalizability. 
These should be regarded as sufficiently restrictive conditions to keep the system's evolution 
in phase space. The Hamiltonian puts more stringent conditions on the admissible class of 
observables, we believe its domain should be dense inside the phase space. Correct normal- 
ization should be then found using the Hamiltonian as a quadratic form on the space of 
these variables and demanding this form to be finite for all physical states. To state our 
conditions we recall that the one-particle Hilbert spaces of bosons are divided into positive 
and negative energy subspaces according to sgn(p_). We assume that the operators which 
act between positive and negative energy subspaces will be Hilbert-Schmidt, whereas the 
operators acting between the same subspaces will be trace class . We write explicitly the 
one for C: C{u-,ui;v-,vi) is trace class for U-V- > 0, and Hilbert-Schmidt for U-V- < 0, 
and the same for the other variables. Note that these are consistent with the constraints on 
the system, which defines the geometry of the phase space. The constraints and the conver- 
gence conditions can be cast into a coherent geometric picture: it defines a a homogeneous 
manifold of Ui{{T-Ca © 'Hb)+, {Ha ® Hb)-), but we will not use it in this work. 

We can rewrite the large- A^j Hamiltonian in terms of these variables, 

H = \j[dp]^^j±^K{p,p) + \j[dp]^^^ 

+ A f[dpdqdsdt]^-^Zl±IzS\Na{p,q)Na{s,t) + Nb{p,q)Nb{s,t)-2C{p,q)C{s,t) 
64^ J\p_q_s_t_\ ^ 



This defines our large- iVj approximation, and in principle we can calculate the equations 
of motion of the basic observables by, 

dO(u,v) , ,_, 

-^-^ = {Oiu,v),H}, (7) 

where O refers to any one of Na, Nb, C, C. The resulting equations are nonlinear integral 
equations and we also have the constraint to satisfy. It would be interesting to study this 
system using a variational ansatz. We will leave the analysis of the full system to a future 
work and look at a linearized version. 



3 Linearization and a possible bound state 

To get a better feeling about the system we can start with a linear approximation. This 
means we should linearize the constraint as well as the equations of motion. The linearization 
of the constraint gives us, for our basic variables, 

(1 + sgn(M_)sgn(f_))A^a(M_,Mi;f_,t;i) = (1 + sgn(M_)sgn(t;_))A^b(M_, Mi; f_, f i) = 
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(sgn(M_) + sgn{v_))C{u^,Ui,v^,Vi) 



0. 



We will be using the last one in our computations, it says that the light-cone momenta should 
be opposite to each other: C{u-,Ui;V-,Vi) = if M_t>_ > 0, and nonzero if u^v^ < 0. The 
same conditions hold for Na and Nf, as well. In the linear approximation we will search 
for a possible bound state of a and b particles. In principle we can compute the linearized 
equations for all the other variables, but they will not lead to a solution for the bound state 
nor a solution for a resonance: they only have scattering states. Therefore we work only with 
the composite C{u,v), let us choose > 0,f_ < 0. The equations of motion of C{u,v) for 
M_ > 0, f _ < 0, in the linear approximation become. 



dC{u,v) 



{C{u,v),H} 



X 



— i- 



[dpdq] 



ml + vl 



q- + v_ 



C{u,v) 

u-\^[pi -qi + Vl - Ui]C{u, v) 



Let us make an ansatz, in the light-front direction we make a 't Hooft like choice with 
respect to the relative momentum variable C, = U-/{u- — f-). This variable now satisfies 
< ( < 1, and we set C(m_, ui; f _, fi) = /(C; tti, f i)e*-^+^ . Notice that P_ = m_ — f_ > 
0, and we introduce a relativistically invariant mass variable fi^ = 2P+|m_ — \ — {ui + 
{—Vi))^, which will be the mass of the bound state (recall that the momentum vi denotes 
an antiparticle with momentum — fi, thus Ui + (— fi) is the total transversal momentum of 
the two particle state). After some manipulations, similar to the ones in 
us an eigenvalue equation for the invariant mass: 



11], this gives 



K + (Ml - Ciui + {-Vi))f ^ ml + {vi - (1 - C)(mi + {-vi))f 



c 

— dr] r i^P^^i^]f(v;pu(ii) 

Stt Jo "^J-oo 



1-C 

<^bi - 9i - (mi + {-Vl))]. 



f{C;ui,Vi) 



(Note that this form reduces to /i = -|- if we set A = 0, and choose the function 
/ properly). We may equivalently use a new set of variables, R = ui — ({ui + (— fi)) 
and Ml + (—Ml), relative transversal light-front momentum and transversal total momentum, 
respectively, instead of the above variables. If we write everything in terms of these new set 
of variables we get. 



/iV~(C ;P,Mi + (-Mi)) 



8n L 



ml + R^ 



c 

[dQMpi 



+ 



m'l 



qi)]fiv;Q,Pi + 



f{C,R,ui + 

-Qi)) 



r^{l-r^)C{l-C) 



Vl)) 

6[pi -qi- (mi + (-Ml))]. 



(The Jacobian of this transformation in the integral is one) 
can be done due to the delta function and we end up with. 



The total momentum integral 



^^f{C-R,ui + {-vi)) 



ml + R^ ^ml + {-Rf 



c 



1-C 



/(C;P,Mi + (-Mi)) 
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[dQ]f{rJ■,Q,u^ + {-v^)) 
^r/(l-r/)C(l-C) 



We see that the total momentum ui + {—Vi) is conserved, thus we can factor out the transver- 
sal center of mass motion by assuming f{(; ui, vi) = /(C, R)g{ui + {—vi)): 



C 



+ 



1-C 



fiCR)- 



X 



167r2 



/oo /•! 
dQ / dr] 
-oo Jo 



(9) 



We can reduce this again to a functional equation for the unknown eigenvalue, by using the 
standart techniques. 



A 



/oo rl 
-oo Jo 



dijdQ 



R^ + ml + jj^rf + [ml — ml — iJ?)r] 



1. 



(10) 



The integrand will have no poles if the quadratic expression involving r] has no real roots, or 
if it has a double root. This is the case if \ma — m},\ < < + m^. We assume this case, 
the last one is clear it says that the bound state cannot be of bigger mass, the other one says 
that the fundamental quanta should be stable against decay, if for example ma > m^ then 
it would be favorable to have ma i— > m(, + /x. Then we can evaluate the integral in any way 
we want, first we take the Q integral, this gives us. 



Jo 



drj 



I li^rf + {ml — ml — ij?)ri + ml 
The next integral can be done and simplified into 



167r 



(11) 



In 



ma + mb + /J. 



ma + mb- n 



IGvr 



(12) 



which is valid when \ma — m},\ < /J' < ma + m^- We may study the small coupling limit of 
this expression. In this case we expect that the bound state mass becomes very close to the 
two mass treshold, then we can write 

+ mfc - /X ^ g-i6,r(m„+m,)/A^ // f« (m„ + m^Ml - 2e-i6'^("*"+"^'')/^), (13) 
ma + mb + ^Ji 

which is consistent if we take X/{ma + m^) << 1. the other extreme is interesting as well, 
fi ^ ma — mb (assuming a is the heavier particle), this implies a critical coupling Ac beyond 
which our methods break down, due to the appearance of a tachyon. 



ma - mb 



In 



mn 



rrib 



or A. 



IGtt 



{ma - mb). 



(14) 



This critical coupling is pushed to higher and higher values if b becomes lighter and lighter 
with respect to the a particle. For any given value of the coupling constant in the interval 
(0, Ac], there is a solution for the bound state energy. Hence we see that there is a composite 
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bound state for these values of the couphng constants. It is not clear what happens beyond 
this value, it is possible that the linear approximation breaks down, it is also possible that 
the large- A^j limit is not a good approximation beyond a certain value. The other possibility 
is that the naive vacuum is not a true vacuum of the quantum theory and we should redefine 
the vacuum of the system. We are not able to analyze these possiblities at the moment. 

Let us compare this with the results we would have found if we looked at a 1 + 1- 
dimcnsional version of the same model. Then the bound state equation could be written 
in terms of the fractional light-cone momentum ( only. There is no transversal component 
and we have only one integral to compute. As a result we find the equation that should be 
satisfied by the eigenvalue /i, (with the condition jm^ — m(,| < /i < -|- rrib), 

Stt 



arctan ( ^ =j + arctan ° =jj = — , (15) 



Amlnil — z^ J Amlml — z^ J ^m\vn^ — z^ 



where = + + z and we wrote for the coupling constant since it has dimensions 
mass-squared. To analyze the behaviour it is more natural to define the dimensionless 
variables, z — z/2mam},, and a — ma/mt,, and rescale the coupling A^ — X^/mamb, 



1 



arctan 



/ o- + z \ , /'l/a + z\-\ 16% 



Note that now z satisfies —l<z<l. If we take the limit z — > — 1+, this corresponds to 
\ma — mfel"*" and the other limit S — > 1~ corresponds to ^ {nia + mb)~ ■ If we assume 
f f=:i 1~ we see that the bound states satisfy a relation 

A^ r A^ 1 
5 ~ 1 - T7^^ or (m„ + mb) 1 - — . (17) 

and if we take A sufficiently small this is consistent. Notice that in 2 + 1 dimensions we 
have an exponential behaviour in the inverse coupling, which is nonanalytic in the coupling 
constant (around zero), as opposed to this power law change. If we look at the other extreme 
we see that there is a finite limit for 5 ^ — 1+, in fact it is equal to 1. This implies a critical 
coupling again, beyond which our methods predict a tachyonic state, 

A^ = IGtt, or A^ = IdTrmamb-- (18) 

This is to be compared with the result in equation (14), which is sensitive to the mass 
difference. 

Let us comment on the convergence conditions in this context. Since we are looking for 
a normalizable solution it looks natural to demand 

rl poo 

/ dC [dR]\fiC,R)\'<oc. (19) 

Jo J-oo 

In fact this is right, and we could see this from our Hilbert-Schmidt condition, 

[du-dv-][duidvi\\C{u-,ui;v-,vi)f < oo. (20) 



u-V-<0 
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If we now make the above change of variables by caUing — — P_ we have 

- f'dC TP-idP-] [[dim + {-v,))dR]\C{P_,C,R,{ui + i-v^m' < oo, (21) 

TT Jo J0+ J 

In our case we are restricting P_ to the surface 2P_P+ = /x^ + {ui + {—Vi)y (for fixed /i, P+), 
this means we should reinterpret the above normalization as 

£ dC J [dR]\fiC,R)\'J ^f-ti!|±i:!^|^(«^ + i-v.mdiu, + i-v,))] < oo, (22) 

(notice that P+ is not allowed to be zero), which means two separate conditions, 

£ dC I [dR]\f{C P)r < oo / ^^^-^^^^^^^^\9{u, + {-v,md{u, + i-v,))] < oo. 

(23) 

The second one simply is a Sobolev type condition, which states that the energy of the 
transversal center of momentum component should also be finite. We see that our solution 
actually satisfies a stronger condition for equations of motion to make sense. 
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